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ABSTRACT. We discuss cosmological inference from galaxy surveys, the X- 
Ray Background (XRB) and the Cosmic Microwave Background (CMB). We 
assume a family of Cold Dark Matter (CDM) models in a spatially flat uni- 
verse with an initially scale-invariant spectrum and a cosmological constant. 
Joint analysis of the CMB and IRAS rcsdhift surveys yields as optimal pa- 
rameters (for fixed U^h 2 = 0.024): 

O m = 1 - Q A = 0.4, h = 0.53, Q = 17 fiK, and fe iras = 1.2. For the above 
parameters the normalisation and shape of the mass power-spectrum are erg = 
0.7 and T = 0.15, and the age of the Universe is 16.5 Gyr. Assuming a 
standard CDM model, the joint IRAS and XRB analysis gives bias parameters 
of &iras = 1-1 and b x (z = 0) = 2.6. When combining CMB and IRAS with 
XRB data we show that standard CDM cannot fit all three data sets. We also 
use these data sets to verify the Cosmological Principle and to constrain the 
fractal dimension of the universe on large scales. 



1 Introduction 

Observations of large scale struc- 
ture (LSS) and the Cosmic Mi- 
crowave Background (CMB) each 
place separate constraints on the val- 
ues of cosmological parameters. 

Estimates derived separately from 
each of these two data sets have 
problems with parameter degener- 
acy. In the analysis of LSS data, 
there is uncertainty as to how well 
the observed light distribution traces 
the underlying mass distribution. 
The light-to-mass linear bias, b, in- 
troduced to account for this un- 
certainty, affects the value of many 
central cosmological parameters, and 
makes any identified optimum de- 
generate. Similarly, on the basis 
of CMB data alone, there is con- 
siderable degeneracy between h = 
J/o/100 km s- 1 Mpc _1 and the en- 
ergy density C2a due to the cosmo- 
logical constant. This leads to poor 
estimation of the baryon (Cb) arl d 
total mass (O m ) densities. 

Several authors (e.g. Gawiser & 
Silk 1998; Eisenstein, Hu & Tegmark 
1998) have recently discussed joint 
analysis of cosmological probes. Here 
we summarize results from Webster 
et al. (1998) and Bridle et al. (in 
preparation) which combine CMB, 



XRB and IRAS data. We present a 
self-consistent formulation of CMB 
and LSS parameter estimation. In 
particular, our method expresses the 
effects of the underlying mass dis- 
tribution on both the CMB poten- 
tial fluctuations and the IRAS red- 
shift distortion. The clustering of 
galaxies in redshift-space is system- 
atically different from that in real- 
space. The mapping between the two 
is a function of the underlying mass 
distribution, in which the galaxies 
are not only tracers, but also veloc- 
ity test particles. The joint analy- 
sis breaks the degeneracy inherent 
in an isolated analysis of either data 
set, and places tight constraints on 
several cosmological parameters. For 
simplicity, we restrict our attention 
to inflationary, Cold Dark Matter 
(CDM) models, assuming a flat uni- 
verse with linear, scale-independent 
biasing. 

2 The CMB 

The compilation of CMB data set 
is described in Hancock et al. (1998) 
and in Webster et al. (1998) and 
is shown in Figure 1. As in Han- 
cock et al. (1998) we form a chi- 
squared x 2 (<^cmb) statistic between 
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Figure 1. A compilation of CMB measurements and the 'best fit' universe 
based on the joint CMB+IRAS estimation, for parameters given in Table 1. 



the flat bandpowers AT; (20 mea- 
surements) for a given set of param- 
eter values (<S cm b)- Since the CMB 
data points were chosen such that no 
two bandpower estimates come from 
experiments which observed overlap- 
ping patches of sky and had over- 
lapping window functions, we may 
consider them as independent esti- 
mates of the CMB power spectrum. 
As the cosmic variance has already 
been taken into account in deriv- 
ing the flat bandpower estimates, the 
likelihood function is given simply by 

^cmb OC e ^ / . 

We assume that the Universe is 
spatially flat, and that there are 
no tensor contributions to the CMB 
power spectrum. We take the pri- 
mordial scalar perturbations to be 
described by the Harrison-Zcl'dovich 
power spectrum for which n B = 1, 
and further assume that the optical 
depth to the last scattering surface is 
zero. 

The normalisation of the CMB 
power spectrum is determined by 
Q, which gives the strength of the 
quadrupolc in fiK. Hereafter fi c d m 
and fit, denote the density of the 
Universe in CDM and the baryons 
respectively, each in units of the crit- 
ical density. Given that we assume a 
flat universe, but investigate models 



where fi m = (fi cdm + fib) < 1 , the 
shortfall is made up through a non- 
zero cosmological constant A such 
that fi A = l-fi m = A/(3H ( 2) . Fur- 
thermore, we restrict our attention to 
models that satisfy the nucleosynthe- 
sis constraint il b h 2 = 0.024 (Tytler 
et al. 1996). Thus we consider the re- 
duced set of CMB parameters 

"cmb = {Q,h, ficdm} ■ (1) 

3 IRAS 

We use the 1.2 Jy IRAS survey 
(Fisher et al. 1995), consisting of 
5313 galaxies, covering 87.6% of the 
sky. 

Here we assume linear, scale- 
independent biasing, where &i ras 
measures the ratio between fluctu- 
ations in the IRAS galaxy distribu- 
tion and the underlying mass density 
field: 

(<Wp)ira S = biras (<Wp) m . (2) 

We note that biasing may be non- 
linear, stochastic, non local, scale de- 
pendent, epoch dependent and type 
dependent (e.g. Dckcl & Lahav 1998, 
Tegmark & Peebles 1998, Blanton et 
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al. 1998, and references therein). For 
a linear bias parameter, b; ras , the ve- 
locity and density fields in linear the- 
ory are linked by a proportionality 
factor /3 iras = r2 m 0i6 /fe iras . 

Statistically, the fluctuations in 
the real-space galaxy distribution 
can be described by a power spec- 
trum, P(k), which is determined 
by the rms variance in the ob- 
served galaxy field, measured for an 
8 ft _1 Mpc radius sphere (erg i ras ) and 
a shape parameter (e.g. T below). 
The observed eg iras is related to the 
underlying erg for mass through the 
bias parameter, such that erg ; ras = 

&iras C8- 

We follow the spherical harmonic 
approach of Fisher, Scharf & Lahav 
(1994). The likelihood of the survey 
harmonics can be calculated as 

Aras °C |A|-3 exp(-5[ 3T S l) .(3) 

Here a is the vector of observed har- 
monics for different radial shells A is 
the corresponding covariance matrix, 
which depends on the predicted har- 
monics (including shot noise). Note 
that the argument of the exponent 
in equation His simply (— x 2 /2), and 
that here the normalisation of the 
likelihood function does depend on 
the free parameters (unlike in the 
CMB likelihood function). 

Since our analysis is valid only in 
the linear regime, we restrict the like- 
lihood computation to Z ma x = 10 
(corresponding to 120 degrees of free- 
dom). The IRAS window functions 
are sensitive to scales k ~ 0.01 — 0.1. 
To summarize, the IRAS likelihood 
function has a parameter vector 

<3iras = {ftras, Cgjras, T} . (4) 

4 Joint analysis CMB+IRAS 

Given the large number of param- 
eters available between the two mod- 
els, it is important both to find links 
for joint optimisation, and to de- 
cide which parameters can be frozen. 
From section 3, we have six vari- 
ables between the two models: {Q, 
h, O cdm , ft ras , cr 8>iras , T}. These can 
be reduced further by expression in 
terms of core cosmological param- 
eters. The IRAS normalisation can 
be calculated as erg = f(Q m ,Q,F), 



while the CDM shape parameter T 
depends on O m h and Q D (Sugiyama 
1995). On the other hand, O m = 

^cdm + ^b; Aras — ^ni 6 /^iras) an d 

f8,iras = o"8^iras- Hence, the final, 
joint parameter space is 

"joint = {h,Q,Qm,b ir!is } ■ (5) 

As the IRAS and CMB probe very 
different scales and hence are as- 
sumed to be uncorrelated, the joint 
likelihood is given by 

ln(£j oint ) = ln(£ cmb )+ln(£ iras ) .(6) 

The joint likelihood (equation ^) 
was maximized with respect to the 
4 free parameters (equation and 
the best fit parameters are shown in 
Table 1. 

Table 1. — Parameter val- 
ues at the joint optimum. The 

68% confidence limits are shown, 
calculated for each parameter by 
marginalising the likelihood over the 
other variables. 



fim 0.39 

h 0.53 

Q (jiK) 16.95 

6iras 1.21 



0.29 < n m < 0.53 
0.39 < h < 0.58 
15.34 < Q < 17.60 
0.98<6 iras < 1.56 



It follows that Q b = 0.085, erg = 
0.67, er 8 , iras = 0.81, T = 0.15, ft ras = 
0.47 and the age of the Universe is 
16.5 Gyr. For this set of parameters, 
we find the values of the reduced x 2 
for the IRAS and CMB data respec- 
tively to be 1.18 and 1.03, confirming 
that both data-sets agree well with 
the models used. Taking the CMB 
and IRAS data together the total re- 
duced x 2 is 1.16. 

We note that our joint IRAS & 
CMB optimal values for erg iras , /9i ras 
and r (Table 1) are in perfect 
agreement with the values derived 
from IRAS alone (Fisher et al. 
1994, Fisher 1994). However at fixed 
f8,iras = 0.69 based on the IRAS cor- 
relation function Fisher et al. (1994) 
found a higher ft ras = 0.94 ± 0.17 
and T = 0.17 ±0.05(1 - cr). 

To obtain 68 per cent confidence 
limits on each of the free parameters 
it is necessary to marginalise over 
the remaining free parameters. The 
marginalised distribution for each 
parameter is shown in Figure 0, 
in which the dashed vertical lines 
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denote the 68% confidence limits 
quoted in Table 1. 

In addition we evaluated the co- 
variance matrix at the joint opti- 
mum. The most strongly correlated 
parameters arc S7 m and h, with 
(normalized) correlation coefficient 
of (-0.82) . 

5 The XRB 

Although discovered in 1962, the 
origin of the X-ray Background 
(XRB) is still unknown, but is likely 
to be due to sources at high redshift 
(for a review sec Fabian & Barcons 
1992). Here we shall not attempt 
to speculate on the nature of the 
XRB sources. Instead, we utilise the 
XRB as a probe of the density fluc- 
tuations at high redshift. The XRB 
sources are probably located at red- 
shift z < 5, making them convenient 
tracers of the mass distribution on 
scales intermediate between those in 
the CMB as probed by COBE, and 
those probed by optical and IRAS 
redshift surveys (see Figure 3). 

The interpretation of the results 
depends somewhat on the nature of 
the X-ray sources and their evolu- 
tion. The rms dipole and higher mo- 
ments of spherical harmonics can be 
predicted (Lahav, Piran & Treyer 
1997) in the framework of growth of 
structure by gravitational instability 
from initial density fluctuations. By 
comparing the predicted multipolcs 
to those observed by HEAOl (Treyer 
et al. 1998) we estimate the ampli- 
tude of fluctuations for an assumed 
shape of the density fluctuations (e.g. 
CDM models). Figure 3 shows the 
amplitude of fluctuations derived at 
the effective scale A ~ 600/i Mpc 
probed by the XRB. The observed 
fluctuations in the XRB are roughly 
as expected from interpolating be- 
tween the local galaxy surveys and 
the COBE CMB experiment. The 

rms fluctuations — on a scale of 
p 

~ 600/i _1 Mpc are less than 0.2 %. 

5.1 Joint IRAS+XRB estima- 
tion 

For simplicity we restricted the 
analysis to the standard CDM model 
(f2 m = l,h = 0.5) with normaliza- 
tion ag = 0.7. Here we assumed a re- 
vised n h h 2 = 0.019 (Buries & Tytlcr 



1998). As the sources of the XRB 
cover a wide range in redshift we as- 
sumed an epoch-dependent biasing 
of the form (Fry 1996): 

b x (z) = b x (0) + z[b x (0) - 1] (7) 

and that the X-ray emissivity varies 
like (1 + z) 2 - 6 out to z max = 6.4 
(Treyer et al. 1998). In the likelihood 
analysis we used the HEAOl data 
with a Galactic mask, and harmon- 
ics 1 < I < 10. We then maximized 
the joint likelihood IRAS+XRB (as- 
sumed to be uncorrelated) with re- 
spect to 2 free biasing parameters, 
and found at the fei ras = 1.1 and 
b x {0) = 2.6. The goodness-of-fit is 
Xl as = 1-22 and X \ rh = 1.01. 

5.2 Joint IRAS+XRB+CMB es- 
timation 

We then added the CMB and 
solved for 3 free parameters: fei ras = 
0.7, M0) = 1.8 and cr 8 = 1.0. While 
the x 2 f° r IRAS and XRB are ac- 
ceptable, the CMB fit is very poor, 
X 2 cmh = 2.7. Hence standard CDM 
cannot fit these 3 data sets simulta- 
neously. We intend to generalise this 
analysis for other models. 

6 Comparison with other 
studies 

The results of the CMB+IRAS op- 
timisation are in reasonable agree- 
ment with other current estimates. 
The relatively low value of f2 m ks 
0.4 is close to that found by oth- 
ers (White et al. 1993, Bahcall et 
al. 1997), and is in line with re- 
cent supernovae results (Perlmut- 
ter et al. 1998). However, given 
the assumption of a flat universe, 
it requires a very high cosmolog- 
ical constant (f2^ = 0.6). Gravi- 
tational lensing measurements have 
constrained Q\ < 0.7 (Kochanck 
1996). Our value for the Hubble con- 
stant, h = 0.53, agrees well with sev- 
eral other measurements (Sugiyama 
1995, Lineweaver et al. 1997), but 
falls at the low end of the gener- 
ally accepted range from local mea- 
surements (Frcedman et al. 1998). 
Assuming the nucleosynthesis con- 
straint Qi,h = 0.024 the optimal 
baryon density is found to be flj, = 
0.085. Our value for the combination 
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Figure 2. The one- dimensional marginalised probability distributions for each 
of the four parameters. The vertical dashed lines denote the 68% confidence 
limits. The horizontal plot limits are at the 99% confidence limits. 



cr 8 f7 m ' 6 = 0.38 is lower than the one 
derived from measurements from the 
peculiar velocity field, crgSl m ' 6 »J 
0.8 (Freudling et al. 1998). Our 
values are closer to the combina- 
tion derived from cluster abundance 
o-$n m - 5 » 0.5 (Eke ct al. 1998). Fi- 
nally, for spatially-flat universes the 
time since the Big Bang for the val- 
ues of our n m and h at the joint op- 
timum is 16.5 Gyr. 

On the IRAS side, ft ras = 0.47 
is in agreement with several other 
measurements (Willick et al. 1997), 
although there are other measure- 
ments which place ftras much higher 
(Sigad et al. 1997). Finally, the 
IRAS mass-to-light bias is seen to be 
slightly greater than unity (6i ras = 
1.2), suggesting that the IRAS galax- 
ies (mainly spirals) arc reasonable 
(but not perfect) tracers of the un- 
derlying mass distribution. On the 
other hand, the XRB sources are 
strongly biased relative to the mass 
(6,(0) = 2.6). 

7 Is the FRW Metric Valid on 
Large Scales ? 

The Cosmological Principle was 
first adopted when observational cos- 



mology was in its infancy; it was then 
little more than a conjecture. Ob- 
servations could not then probe to 
significant rcdshifts, the 'dark mat- 
ter' problem was not well-established 
and the CMB and the XRB were still 
unknown. If the Cosmological Prin- 
ciple turned out to be invalid then 
the consequences to our understand- 
ing of cosmology would be dramatic, 
for example the conventional way 
of interpreting the age of the uni- 
verse, its geometry and matter con- 
tent would have to be revised. There- 
fore it is important to revisit this 
underlying assumption in the light 
of new galaxy surveys and measure- 
ments of the XRB and CMB. The 
question of whether the universe is 
isotropic and homogeneous on large 
scales can also be phrased in terms 
of the fractal structure of the uni- 
verse. A fractal is a geometric shape 
that is not homogeneous, yet pre- 
serves the property that each part is 
a reduced-scale version of the whole. 
If the matter in the universe were 
actually distributed like a pure frac- 
tal on all scales then the Cosmolog- 
ical Principle would be invalid, and 
the standard model in trouble. As 
shown in Figure 3 current data al- 
ready strongly constrain any non- 
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Figure 3. A compilation of rms density fluctuations, (^) 2 ~ k 3 P(k), on 
different scales from various observations: a galaxy survey, the X-ray Back- 
ground and Cosmic Microwave Background experiments. The crosses repre- 
sent constraints from the XRB HEAOl quadrupole (Lahav et al. 1997, Treyer 
et al. 1998). The top and bottom crosses are estimates of the amplitude of the 
power- spectrum at fc _1 ~ GOOh^ 1 Mpc, assuming CDM power-spectra with 
shape parameters F = 0.2 and 0.5 respectively, and an Einstein-de Sitter uni- 
verse. The fractional error on the XRB amplitudes ( due to the shot-noise of 
the X-ray sources) is about 30%. The solid and dashed lines correspond to the 
standard CDM power-spectrum (with shape parameter V = 0.5j and a 'low- 
density' CDM power- spectrum (with T = 0.2), respectively, assuming erg = 1 
in both cases. The open squares at small scales are estimates of the power- 
spectrum from 3D inversion of the angular APM galaxy catalogue (Baugh & 
Efstathiou 1994). The elongated box at large scales represent the COBE 4~yr 
CMB measurement. The COBE box corresponds to a quadrupole Q=18.0 fj,K 
for a Harris on- Zeldovich mass power-spectrum, via the Sachs- Wolfe effect, 
or as = 1.4 for a standard CDM model (Gawiser and Silk 1998). 



uniformities in the galaxy distribu- 
tion (as well as the overall mass dis- 
tribution) on scales > 300 ft -1 Mpc. 

If we count, for each galaxy, the 
number of galaxies within a distance 
R from it, and call the average num- 
ber obtained N(< R), then the dis- 
tribution is said to be a fractal of cor- 
relation dimension D2 if N(< R) oc 
R D2 . Of course D2 may be 3, in 
which case the distribution is homo- 
geneous rather than fractal. In the 
pure fractal model this power law 
holds for all scales of R. 

The fractal proponents 

(Pietronero et al. 1997) have esti- 
mated D2 ~ 2 for all scales up to ~ 
1000 h^ 1 Mpc, whereas other groups 
have obtained scale-dependent val- 



ues (for review see Wu, Lahav & Rees 
1998, and references therein). 

If we assume homogeneity on large 
scales, then we have a direct mapping 
between correlation function £(r) (or 
the Power-spectrum) and D2- For 
£(r) oc r ' it follows that D2 = 3 — 7 
if £ > 1, while if £(r) = then 
D2 = 3. We note that it is inap- 
propriate to quote a single crossover 
scale to homogeneity, for the transi- 
tion is gradual. 

Direct estimates of D2 are not 
possible for much larger scales, but 
we can calculate values of D2 at 
the scales probed by the XRB and 
CMB by using CDM models nor- 
malised with the XRB and CMB as 
described above. The resulting val- 
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ues are consistent with D2 = 3 to 
within 10~ 4 from the XRB on scales 
~ 500ft _1 Mpc and 2 X 10~ 5 from 
the CMB on ~ 1000 h~ x Mpc (Wu et 
al. 1998). Isotropy does not imply ho- 
mogeneity, but the near-isotropy of 
the CMB can be combined with the 
Copernican principle that we are not 
in a preferred position. 



8 Discussion 

The near future will see a dra- 
matic increase in LSS data (e.g. 
the PSCZ, SDSS, 2dF and 6dF sur- 
veys) and detailed measurements of 
the CMB fluctuations on sub-degree 
scales (e.g. from the Planck Surveyor 
and MAP satellites). These will al- 
low more accurate parameter esti- 
mation and exploration of a wider 
range of models. In particular we em- 
phasise that the naive linear bias- 
ing should be generalised to more re- 
alistic scenarios. Other cosmological 
probes (e.g. Supernovae, clusters, pe- 
culiar velocities and radio sources) 
can be added to the analysis to set 
tighter constraints on the parameter 
space. 
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